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I.  INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

The  need  for  improved  modeling  of  aircraft  aerodynamic 
characteristics  has  been,  and  continues  to  be,  apparent  in 
numerous  areas  of  great  operational  importance.  Reliable  and 
accurate  flight  simulators  for  pilot  training,  verification  of 
design  methods  for  specification  of  aircraft  characteristics, 
and  the  off-line  development  of  mission  profiles  that  make 
optimum  use  of  the  airplane's  capabilities  (without  exceeding 
its  limitations)  are  three  such  areas.  In  general,  there  is 
also  a  need  for  improved  understanding  of  an  airplane's  intrin¬ 
sic  aerodynamic  characteristics  which,  if  filled,  will  lead  to 
improved  cost  effectiveness  and  operational  safety. 

The  ability  to  process  large  quantities  of  data  from 
flight  tests  has  led  to  the  parallel  development  of  data  proces¬ 
sing  algorithms  which  greatly  increase  the  amount  of  useful 
knowledge  that  can  be  extracted  from  the  data.  These  algorithms 
based  on  dynamical  and  statistical  principles,  yield  very  pre¬ 
cise  information  about  the  characteristics  of  the  data  and  the 
system  that  produced  it. 

1 . 2  BACKGROUND 

The  present  study  represents  the  application  of  an  advanced 
system  identification  methodology  which  has  been  developed  speci¬ 
fically  for  estimating  high-angle-of-attack  aerodynamics  stability 
and  control  coefficients  from  flight  data.  The  development  of 
this  methodology  has  been  documented  extensively  in  previous 
technical  reports  and  technical  journals,  and  is  summarized 
historically  as  follows. 


the  algorithms  on  linear  aircraft  response  data.  This  was  done 
in  a  two-step  process  by  processing  F-14  data  from  the  Naval  Air 
Test  Center  and  subsequently  installing  the  software  at  NATC. 

The  identification  program  installed  at  NATC  had  the  ability  to 
identify  both  process  and  measurement  noise  and  the  operation  of 
these  programs  is  detailed  in  Ref.  5.  The  program  was  further 
applied  to  other  aircraft  by  NATC  personnel  as  discussed  in 
Refs.  6  and  7. 

(4)  Application  of  the  software  to  simulated  and  actual 
flight  test  data  indicated  the  need  for  a  near  real-time  identi¬ 
fication  algorithm  to  provide  rapid  estimates  of  aircraft  linear 
stability  and  control  derivatives.  A  theoretical  basis  for  such 
reduced  computation  time  was  established  with  further  research, 
resulting  in  two  complementary  algorithms.  The  first  algorithm 
was  based  on  the  theory  of  Luenberger  canonical  forms  for  linear 
systems  and  reduced  execution  time  of  the  linear  maximum  likelinood 
algorithm  by  a  factor  of  ten  (depending  on  the  number  of  param¬ 
eters  to  be  identified).  The  second  algorithm  was  based  on  a  new 
formulation  of  an  instrumental- variables  approach  [8]  and  pro¬ 
vided  true  real-time  estimation  capability  for  linear  aircraft 
stability  derivatives.  The  theoretical  development  and  implemen¬ 
tation  of  these  two  algorithms,  with  application  to  simulated  and 
actual  flight  test  data,  is  presented  in  Ref.  5. 

(5)  The  developments  of  the  previous  phases  had  demonstrated 
the  further  requirements  for  expansion  of  the  theoretical  base 
upon  which  nonlinear  system  identification  was  based.  Specific 
research  into  both  the  nonlinear  model  structure  determination 
and  input  design  issues  was  performed,  the  results  of  which  are 
reported  in  Ref.  6. 

The  results  of  this  report  are  the  application  of  the  identi¬ 
fication  methodology  to  actual  nonlinear  aerodynamic  regime  data. 
The  data  for  this  analysis  were  provided  by  the  Naval  Air  Develop¬ 
ment  Center.  The  specific  objectives  are  listed  in  Table  1.1. 


TABLE  1.1.  OBJECTIVES 


•  DETERMINE  NONLINEAR  T-2  AERODYNAMIC 
COEFFICIENTS  AND  MODEL  STRUCTURES  IN 

HIGH  ANGLE  OF  ATTACK  AND  SIDESLIP  REGIMES. 

•  EVALUATE  MODEL  CHARACTERISTICS  AGAINST 
EXISTING  DATA. 

•  EVALUATE  MODEL  PARAMETERS. 

•  DEMONSTRATE  ADVANCED  NONLINEAR  IDENTI¬ 
FICATION  TECHNIQUES. 


1 . 3  SUMMARY 

Subsequent  sections  of  this  report  are  organized  as  follows. 
Section  II  discusses  the  T-2C  flight  data  analyzed  in  this  effort, 
and  describes  wind  tunnel  data  and  theoretical  aerodynamic  tech¬ 
niques  used  as  a  guide  in  model  structure  determination.  Section 
III  describes  the  data  preparation  techniques  employed,  and  the 
application  and  findings  of  the  model  structure  determination 
and  model  parameter  estimation  techniques.  Section  IV  summarizes 
the  study  results  and  presents  conclusions  and  recommendations. 

Various  detailed  technical  data  are  presented  in  the 
appendices. 
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II.  T-2C  AERODYNAMIC  DATA  BASE 


2.1  OVERVIEW 

This  section  discusses  the  flight  test  data  used  in  the 
determination  of  T-2C  high  angle  of  attack  aerodynamic  coefficients 
through  model  development  and  parameter  estimation,  and  data  used 
to  support  modeling  efforts  and  give  additional  understanding  of 
aircraft  characteristics.  Wind  tunnel  data  and  theoretical  aero¬ 
dynamic  calculations  were  used  to  validate  each  other  and  provide 
a  reference  of  comparison  for  the  identified  models.  In  this  way, 
as  thorough  an  understanding  of  T-2C  aerodynamics  as  could  be 
obtained  was  developed  to  guide  the  aerodynamic  model  analysis. 

2.2  FLIGHT  TEST  DATA 

The  flight  test  data  analyzed  in  this  study  were  taken  from 
flights  made  between  5  November  1975  and  26  August  1976.  The 
aircraft  configuration  during  these  flights  is  summarized  in 
Table  2.1;  it  represents  the  normal,  clean  aircraft  configur¬ 
ation.  Detailed  aircraft  geometry  from  Ref.  9  is  given  in 
Appendix  A. 

The  initial  condition  for  all  flight  maneuvers  was  nominally 
as  summarized  in  Table  2.2;  small  deviations  from  the  desired  trim 
condition  were  accounted  for  in  the  data  analysis  procedure. 
Significant  to  the  goal  of  nonlinear  aerodynamic  coefficient 
determination  are  the  high  initial  value  of  and  the  Reynolds 

number.  The  trim  angle  of  attack  is  just  3.5°  below  the  stall 
angle  of  attack  predicted  by  both  wind  tunnel  data  and  theoretical 
methods  (to  be  discussed  below),  and  the  trim  lift  coefficient  of 

ClTRIM  *  *93  is  78  Percent  of  trim  cLmax*  Because  the  flight 
and  wind  tunnel  Reynolds  numbers  are  quite  close  (7.5  x  106  vs 
4  x  IQ**,  respectively),  any  differences  between  flight  and  wind 
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TABLE  2.1. 


T-2  FLIGHT  AIRCRAFT  CONFIGURATION  SUMMARY 


•  CLEAN  CONFIGURATION 

-  FLAPS  NEUTRAL  (3°  UPRIG) 
-AILERONS  NEUTRAL  (3°  UPRIG) 

-  GEAR  UP 

-  SPEED  BRAKE  RETRACTED 

•  NEARLY  CONSTANT  WEIGHT 

-  1 1 ,400  LBS  +  400  LBS 

•  THROTTLES  FIXED  DURING  MANEUVERS 

•  INERTIAS  (SLUG-FT2) 

Ix  =  9000 

I  =  14,600 
Iz  =  19,000 


TABLE  2.2  FLIGHT  TEST  CONDITION 


V  =  300  FT/ SEC  (TAS) 

MACH  =0.3 

h  *  25,000  FT  (PRESSURE) 

W  =  11,400  LB 

yTRIM*  °’° 

REYNOLDS  NUMBER  =  7.5 x  106 

“TRIM  =  10-5° 

CLTRIM  =  '93 

Clmax  ’  1,2 


NOTE:  DATCOM  ANALYSIS  SHOWS  NEGLIGIBLE  COEFFICIENT 
VARIATION  DUE  TO  CHANGE  FROM  WIND  TUNNEL 
RN  =  4  x  10®  TO  FLIGHT  R^  =  7.5  x  106 


tunnel  characteristics  are  most  likely  due  to  effects  from 
other  than  Reynolds  number.  Theoretical  studies,  using  DATCOM 
predictions,  verify  this  conclusion. 


Large  perturbation  maneuvers  were  initiated  at  this  flight 
condition  by  various  deflections  of  one  or  more  control  surfaces, 
with  the  goal  of  forcing  the  aircraft  into  the  stall  and  post 
stall  regions  of  angle  of  attack  with  sizeable  sideslip  excursions 
Seventeen  maneuvers,  shown  in  Table  2.3,  were  selected  for  analysi 
in  this  study.  The  types  of  control  inputs  employed  are  sketched 

TABLE  2.3.  T-2C  FLIGHT  TEST  CONDITIONS  USED  FOR 
NONLINEAR  AERODYNAMIC  IDENTIFICATION 


DATE 

FLT.  TEST 
NO. 

INPUT 

NOV.  5.  1975 

S 

£e  AFT,  FULL  6r 

9 

6e  AFT,  FULL 

DEC.  10.  1975 

4 

&e  RANDOM 

8 

(t  RANDOM 

14 

SEQUENTIAL  DOUBLETS 

DEC.  29,  1975 

2 

6e  DOUBLET 

14 

ir  RANDOM 

17 

SEQUENTIAL  RANDOM  (et«  0°) 

AUG.  26,  1976 

1 

6g  RAMP  +  SINE 

2 

SEQUENTIAL  RANDOM  (£t*  -5°) 

3 

SEQUENTIAL  RANDOM  (6t*+5°) 

4 

ie  AFT  +  PULSES 

5 

LIMIT  CYCLE 

6 

SHALLOW  BANK  TO  STALL 

7 

STEEP  BANK  TO  STALL 

8 

PULL  UP  TO  STALL 

9 

COORDINATED  SPIN  ENTRY 

in  Figure  2.1.  These  inputs  were  made  by  the  pilot  and  are  there¬ 
fore  neither  repeatable  nor  precise  with  regard  to  their  spectral 
characteristics,  but  are  generally  effective  in  creating  large- 
amplitude  motions.  A  useful  indication  of  excitation  effectiveness 
is  a  raw-data  a- B  plot  (histogram)  such  as  shown  in  Figure  2.2, 
which  provides  knowledge  about  the  type  of  information  that  can 
be  obtained  from  the  maneuver.  Histograms  of  other  flight 
maneuvers  are  shown  in  Section  III.  Overall,  the  maneuvers 
encompassed  angles-of-attack  from  -4°  to  +30°  and  angles-of- 
sideslip  in  excess  of  +^20°. 


The  apparent  effect  of  stall  entry  dynamics  (i.e.,  deceler¬ 
ation  rate  and  pitch  rate)  on  Cl^ax  is  not  an  issue  since  the 
identification  of  (see  Section  III)  is  based  on  accelerometer 
and  airspeed  measurements.  With  this  approach  the  problem  of 
relating  to  air  speed  via  the  weight  of  the  aircraft  for  non 
"lg”  flight  is  avoided.  In  general,  because  the  aerodynamic 
parameters  are  identified  through  the  minimization  of  accelero¬ 
meter,  airdata  and  rate  gyro  measurements,  it  is  extremely 
desirable  for  the  system  identification  techniques  described 
in  this  report  to  maximize  the  dynamic  nature  of  the  flight  test 
maneuver . 


Flight  data  were  simultaneously  recorded  on  board  and  tele¬ 
metered  to  a  ground  receiving  station.  Flight  instrumentation 
provided  6-DOF  measurement  capability  for  both  rates  and 
acceleration.  The  types  and  locations  of  the  sensors  are  shown 
in  Table  2.4,  and  their  respective  ranges  shown  in  Table  2.5. 
Note  that  body-axis  angular  accelerations  are  measured  by  using 
translational  accelerometers  at  the  airplane's  nose,  tail,  and 
wingtips . 


2.3  WIND  TUNNEL  DATA 

An  extensive  wind  tunnel  investigation  of  the  high-a  char¬ 
acteristics  of  the  T-2C  is  reported  in  Ref.  9,  presenting  esti¬ 
mates  of  static  aerodynamic  coefficients  and  control  effectiveness 
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1.  J4  RANDOM 


2.  5#  RANDOM 

3.  5r  RANDOM 

4.  SEQ.  RANDOM 

5.  SEQ.  DOUBLET 


4jy\r 
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hr- 
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n 
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Figure  2.1.  Control  Inputs  for  Parameter  Identification 


NOTES: 


1)  5,  4fT.  fUU  !,  NRNEUNER 

2)  NIMERS  ON  PLOT  SHOW  THE  «N(|  OF  ORTA  POINTS 
rOH  INDICATED  i/a  REGION 


Figure  2.2.  Example  a/0  Histogram 
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TABLE  2.4.  INSTRUMENT  LOCATIONS  ON  YT-2B  FOR  AIRFRAME 
DYNAMICS  IDENTIFICATION  PROGRAM 


FUNCTION 

NAME 

FUSELAGE 

STATION 

WING 

STATION 

WATER 

LINE 

nx 

GREENLEAF  SERIAL 
1655000- A  #674 

37. 

-4. 

-11. 

nz  (c.g. ) 

GIANNINI  #1816 

241 17P-3. 5-20 

207. 

-1. 

-21. 

n2  NOSE  (FOR  q) 

GIANNINI  #1819 
24117P-3.5-20 

34. 

-6. 

-11. 

nz  TAIL  (FOR  q) 

GIANNINI  #1837 
24117P-3.5-20 

366.5 

0 

+9. 

q 

STAHAM  ANG.  ACC. 
AA196-8-350 

230. 

0 

-23.5 

e 

GIANNINI  MODEL 

SERIAL  5812  3416 

56. 

0 

5. 

GIANNINI  MODEL 

SERIAL  3812  3416 

56. 

0 

5. 

r 

NORDEN  RG  228 

SERIAL  #185 

47. 

0 

3.5 

P 

HUMPHREY  G20-1021-00 
SERIAL  #828 

52. 

0 

2.5 

q 

NORDEN  RG  228 

SERIAL  #183 

196. 

0 

-22.5 

ny  NOSE  (FOR  r) 

EDELCLIFF  MODEL  7-30 
SERIAL  #366.5 

31.5 

-6. 

-11. 

ny  (c.g.) 

EDELCLIFF  MODEL  7-30 
SERIAL  #4444 

204. 

-1.5 

-19. 

nz  WING  TIP  (FOR  p) 

GIANNINI  #1685 

241 17P-3. 5—20 

249. 

206.5 

15. 

n2  WING  TIP  (FOR  p) 

GIANNINI  #1681 

241 1 7P-3. 5-20 

249. 

-206.5 

15. 

a 

VANE 

-62.3 

0 

-21 

B 

VANE 

-65.9 

0 

-21 

NOTE:  NEGATIVE  WING 
POSITIONS  ARE 

STATIONS  ARE  TOWARD  LEFT  WING.  NEGATIVE  WATER  LINE 

BELOW  FUSELAGE  REFERENCE  LINE. 
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TABLE  2.5.  MEASUREMENT  RANGES  FOR  YT-2B/C  TRANSDUCERS 


MEASUREMENT 

RANGE 

n 

z 

±3.5g 

"x 

±1 .0g 

ny 

±1 .0g 

q 

±20  DEG/SEC 

e 

±15  DEG 

r 

±20  DEG/SEC 

q 

±450  DEG/SEC2 

* 

±90  DEG 

p 

±45  DEG/SEC 

V 

0  -  500  KTS  (COARSE) 

0-50  ICTS  (FINE  RESETTABLE) 

a 

±10  DEG 

ABOUT  SET  POINT 

0 

±10  DEG 

parameters.  The  overall  ranges  of  simultaneous  a  and  6 
settings  in  the  reference  study  were  -8°  to  +83°  and  -10°  to  +30 
respectively,  accomplished  by  the  use  of  three  different  model 
support  arrangements,  each  allowing  examination  of  a  particular 
range  of  a  within  which  sting  interference  with  the  model  was 
small . 


Figures  2.3  and  2.5  (Ref.  9)  show  basic  lift,  drag,  and  pitching 
moment  coefficient  data  for  8=  0  over  the  entire  a  range. 

Figures  2.6  and  2.7  show  the  variation  of  rolling  and  yawing  moment 
coefficients  with  sideslip  angle  at  various  angles  of  attack. 

These  wind  tunnel  data  were  very  useful  in  the  present 
investigation.  Through  study  of  the  nonlinear  aerodynamic 
characteristics  revealed  in  the  data,  the  nonlinear  functional 
forms  of  the  coefficients  were  obtained  with  less  effort  during 
the  model  structure  determination  stage.  Knowledge  of  key  per¬ 
formance  and  stability  parameters,  such  as  Cr.,.v,  ar_  ,  and  C 

MAX’  lLmax’  ma 

allowed  faster  interpretation  of  flight  data  characteristics  and 
permitted  reasonableness  checks  on  the  hypothesized  functional 

forms  of  the  coefficients.  Any  uncertainty  in  the  wind  tunnel 
data  did  not  substantially  impair  its  usefulness  in  this  effort, 
as  it  was  demonstrated  in  the  work  to  follow  that  the  MSD  and 
parameter  estimation  techniques  require  only  basic,  qualitative 
a  priori  knowledge  of  aircraft  characteristics. 

2.4  THEORETICAL  AERODYNAMIC  DATA 

Theoretical  methods  for  computing  aircraft  aerodynamic 
characteristics  are  necessary  in  order  to  confirm  wind  tunnel 
data  and  provide  estimates  of  aerodynamic  parameters  not  amenable 
to  experimental  determination.  It  is  desirable  that  the  capa¬ 
bilities  of  theoretical  methods  be  well  understood,  as  they  provide 
an  opportunity  to  reduce  the  extent  of  wind  tunnel  testing  required 
to  gain  an  understanding  of  basic  configuration  characteristics. 
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ANGLE  Of  ATTACK,  a,  DEG 


Figure  2.3.  T-2C  Lift  Coefficient  Versus  Angle  of  Attack 
(Wind  Tunnel  Data  from  NADC-73259-30) 


note: 

t  8*0° 


*  0* 
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Figure  2.4. 


T-2C  Drag  Coefficient  Versus  Angle  of  Attack 
(Wind  Tunnel  Data  from  NADC- 73259- 30) 
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Theoretical  aerodynamic  methods  exist  in  a  wide  variety  of 
forms  extending  from  the  implementation  of  closed-form,  configuration- 
dependent  equations  to  elaborate  potential  flow  methods  using  vortex 
lattices  or  source  panels  to  represent  the  aircraft.  The  latter 
approach  is  typically  very  complex  and  costly  in  terms  of  computer 
time,  while  the  former  is  usually  unable  to  handle  the  critical 
matter  of  downwash  distribution  and  the  contributions  of  horizontal 
and  vertical  tail  loads.  A  compromise  approach  has  been  reached  in 
the  methods  of  the  DATCOM,  which  is  based  on  wind  tunnel 
data  unified  by  proven  theoretical  results.  This  method  gives 
aerodynamic  characteristics  of  a  wide  range  of  possible  configura¬ 
tions  to  the  accuracy  required  for  preliminary  design  studies.  The 
utility  of  the  procedure,  which  is  thorough  but  laborious,  has  been 
enhanced  manyfold  by  its  incorporation  in  a  digital  computer  pro¬ 
gram.  The  user,  by  inputting  only  basic  configuration  data,  may 
thus  easily  obtain  estimates  of  static  and  dynamic  aerodynamic 
coefficients  and  derivatives,  several  extending  into  the  nonlinear 
range  of  a  and  8.  The  burden  is  on  the  user,  however,  to 
understand  the  limitations  of  the  methods  and  recognize  those 
areas  in  which  the  predictions  may  be  inaccurate. 

One  of  the  key  benefits  of  this  method  is  the  ability  it 
brings  to  examine  the  effects  of  small  changes  from  the  wind  tunnel 
model  configuration.  The  effects  of  changing  Reynolds  number, 
for  instance,  or  of  removing  the  3°  uprig  in  the  T-2C  flaps  and 
ailerons,  were  studied  in  this  effort. 

For  a  relatively  simple  aircraft  configuration  such  as  the  T-2C, 
the  DATCOM  predictions  for  basic  characteristics  possess  a  high 
degree  of  credibility  at  angles-of-attack  below  stall.  Figure  2.8 
shows  the  lift  curve  prediction  compared  to  wind  tunnel  data.  Lift 
curve  slope,  stall  angle  of  attack,  and  post-stall  lift  behavior 
are  all  in  good  agreement  with  wind  tunnel  findings .  The  good 
prediction  of  stall  a  is  valuable  because  of  its  importance  in 
spline  curve  fitting  for  model  structure  determination.  Figure  2.9 
shows  the  pitching  moment  curve,  which  displays  good  agreement  with 
wind  tunnel  data  below  stall.  Among  the  other  results  obtained 
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Figure  2.8.  T-2C  Comparison  of  Datcom  and  Wind 
Tunnel  Pitching  Moment  Prediction 


Figure  2.9.  T-2C  Comparison  of  natcom  and  Wind 
Tunnel  Lift  Prediction 


are  lateral -directional  damping  and  coupling  derivatives,  of  which 
four  important  coefficients  are  shown  in  Figures  2.10  to  2.13.  The 
DATCOM  predictions  are  compared,  in  these  figures,  to  theoretical 
estimates  computed  by  aerodynamic  strip  theory  methods  as  described 
in  Ref.  10.  The  general  agreement  shown  here  is  encouraging. 

With  the  knowledge  thus  obtained,  from  wind  tunnel  and 
theoretical  methods,  of  aerodynamic  parameters  of  known  importance 
to  flight  dynamics  analysis,  the  evaluation  of  flight  data  is 
begun  to  verify  or  alter  these  parameters  and  search  for  new,  non¬ 
linear  parameters  whose  determination  can  result  in  improved 
dynamic  predictions. 


Figure  2.10.  Comparison  of  CA  Predictions 


Figure  2.  11.  Comparison  of  C  Predictions 
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Figure  2.13.  Comparison  of  Cn  Predictions 
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III.  T-2C  SYSTEM  IDENTIFICATION  RESULTS 


3.1  OVERVIEW 

This  section  discusses  the  application  of  the  system  identifi¬ 
cation  process  to  T-2  flight  test  data  and  presents  the  results  of 
the  nonlinear  aerodynamic  modeling  analyses.  The  objective  of  this 
effort  was  to  obtain,  from  flight  data,  a  set  of  nonlinear  func¬ 
tional  relationships  describing  the  aerodynamic  characteristics  of 
the  airplane,  in  the  form  of  force  and  moment  coefficients  and 
stability  derivatives  expressed  as  functions  of  various  combina¬ 
tions  of  angle  of  attack  and  sideslip,  rotational  rates,  and  control 
deflections.  Dynamic  nonlinear  effects  such  as  aerodynamic 
hysteresis  were  also  sought. 

The  steps  in  the  integrated  system  identification  procedure 
for  the  T-2C  are  summarized  in  Figure  3.1  Five  basic  steps  are 
shown : 

(1)  Data  preparation  (reconstruction  and  filtering] . 

(2)  Model  structure  determination. 

(3)  Model  parameter  estimation. 

(4)  Model  verification. 

(5)  Flight  test  planning. 

The  first  tour  steps  comprise  the  scope  of  the  current  con¬ 
tract  and  are  discussed  in  detail  in  this  chapter.  The  last 
step,  Flight  Test  Planning,  is  reported  in  Ref .  11. 

3.2  FLIGHT  DATA  PROCESSING 

A  significant  element  of  the  system  identification  procedure 
is  the  preprocessing  of  the  flight  test  data.  The  overall 
objective  of  this  step  is  to  develop  a  set  of  kinematically 
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•  STABILITY  ANO  CONTROL  COEFFICIENTS 

•  ERROR  CALIBRATION 


Figure  3.1.  Integrated  System  Identification  Procedure 
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consistent  measurements  for  the  model  structure  determination  and 
parameter  identification  algorithms.  Measurement  consistency  has 
a  significant  impact  on  parameter  identification  accuracy  since 
unaccounted  for  inconsistencies  will  bias  the  parameter  estimates. 
The  steps  followed  in  the  preprocessing  of  T-2C  flight  test  data 
are  now  described. 

3.2.1  Data  Response  Review 

The  time  responses  were  examined  to  establish  the  validity 
of  the  SCI-created  data  files,  to  classify  the  flight  conditions 
in  terms  of  aerodynamic  regimes  which  were  encountered,  and  to 
provide  a  preliminary  assessment  of  data  quality  and  consistency. 
From  the  preliminary  assessment,  requirements  for  filtering  noisy 
data  were  established,  wild  points  and  data  drop-outs  were  visually 
isolated  and  corrected,  and  obvious  measurement  sign  errors  were 
identified. 

Figure  2.2  illustrates  one  of  the  methods  used  to  classify  a 
test  condition.  This  figure  is  a  histogram  of  angle  of  attack 
and  sideslip.  The  numbers  in  the  plot  indicate  the  number  of  data 
points  within  a  particular  a/6  region. 

3.2.2.  Data  Reconstruction 


Data  reconstruction  is  an  essential  subtask  of  the  flight 
data  processing  procedure  and  it  comprises  three  steps. 

First,  since  the  equations  of  motion  which  are  processed  by 
the  identification  algorithms  assume  a  rigid  body  formulation, 
flight  measurements  representative  of  the  motion  at  the  center 
of  gravity  are  required.  Table  3.1  lists  the  equations  which  are 
used  for  the  reconstruction  of  equivalent  c.g.  measurements. 

Second,  noisy  data  are  filtered  to  provide  numerical  condi¬ 
tioning.  The  T-2C  flight  test  data  were  processed  with  a  Martin- 
Graham  digital  filter.  This  digital  filter  has  the  property 
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TABLE  3.1.  EQUATIONS  USED  TO  RECONSTRUCT  C.G.  MEASUREMENTS 


-  VERTICAL  &  RATE  GYROS  DO  NOT  REQUIRE  CORRECTION 

-  AIR  DATA  SYSTEM 

*  AIRSPEED  DOES  NOT  REQUIRE  CORRECTION 
acg  "  Tan  1  CTan  aM  +  V^Cosa,,,  CosB^ 

*  Bcg  *  Sin_1  [Cos  aM  Tan  ®M  ' 

-  ACCELEROMETER 

*  ax  =  a  +  (Q2  +  R2)1  +  (R  -  PQ)£  -  (PR  +  QH 

xcg  xa  ya  za 

*  av  =  a  -  (PQ  +  RU  +  (P2  +  P2U  +  (P  -  RQ)t 

ycg  y  xa  ya  za 

*  a7  =  a  +  (Q-  PRU  -  (RQ  +  P)i  +  (Q2  +  P2)d7 

zcg  z  a  ya  za 


of  not  changing  either  the  magnitude  or  phase  content  of  the  data 
up  to  the  cutoff  frequency  and  then  truncating  the  signal  for 
frequencies  greater  than  the  termination  frequency.  For  the  T-2C 
application,  the  cutoff  frequency  was  set  at  2Hz  and  the  termination 
frequency  at  4Hz.  Since  the  selected  cutoff  frequency  for  the 
Martin-Graham  filter  is  much  greater  than  the  T-2’s  rigid  body 
frequencies  of  motion,  its  application  will  not  change  the  infor¬ 
mation  content  of  the  flight  test  data. 

The  third  step  involves  reconstruction  of  measurements  found  (as 
a  result  of  performing  a  kinematic  consistency  evaluation)  to  be  un- 
useable,  and  by  combining  actual  measurements  to  formulate  needed 
values.  An  example  of  the  first  type  of  reconstruction  would  be 
to  define  roll  acceleration  by  taking  the  time  derivative  of  roll  rate 
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Deriving  roll  acceleration  by  taking  the  difference  between  normal 
acceleration  signals  from  right  and  left  wing  tip  accelerometers 
is  an  example  of  the  second  type  of  reconstruction.  Both  types 
of  reconstruction  were  used  for  processing  the  T-2C  flight  test  data. 

3.2.3  Data  Consistency  Analysis 

The  kinematic  consistency  of  the  sensor  signals  is  evaluated 
as  a  means  of  highlighting  potential  measurement  system  errors. 
Measurement  system  errors  can  result  from  a  number  of  possible 
sources.  Bias  and  scale  factor  errors  could  result  from  using 
an  incorrect  calibration  curve  or  from  the  sensor's  power  supply 
not  operating  at  its  design  condition.  Since  a  body  axis  frame 
of  reference  is  used  for  solving  the  equations  of  motion,  the 
orientation  of  the  sensor  system  relative  to  the  body  axis  system 
must  be  accurately  known.  Instruments,  such  as  a  vertical  gyro 
(VG) ,  can  produce  erroneous  measurements  as  a  result  of  their 
operational  characteristics.  For  example,  if  the  VG ' s  erection 
circuit  is  allowed  to  operate  during  accelerated  flight,  the 
attitude  measurements  will  be  referenced  to  a  false  vertical, 
and  thus  will  be  in  error.  This  type  of  problem  would  be  most 
significant  for  the  identification  of  low  frequency  modes  of 
motion.  Another  source  of  inconsistency  exists  when  inertial 
velocity  components  are  not  measured  (this  is  the  case  for  the 
T-2C  program).  Since  the  equations  of  motion  are  written  in 
terms  of  inertial  velocity  components,  substitution  of  air 
referenced  components  (i.e.,  V^.  cosB  cosa  for  u)  will  produce 

biased  answers  when  winds  are  present  and  not  accounted  for. 

An  example  of  a  data  consistency  check  is  shown  in  Figure  3.2, 
where  the  variation  of  q  with  time  does  not  agree  with  dq/dt 
(q  from  a  rate  gyro)  or  with  q  computed  from  0  and  <J>  from 
vertical  gyro).  The  latter  two  approaches  agree  with  each  other, 
so  the  conclusion  is  that  the  q  measurement  is  faulty  and  should 
be  replaced  with  one  of  the  other  (derived)  measurements. 
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3.2.4  T-2C  Flight  Data  Quality 


As  a  result  of  processing  the  T-2C  flight  test  data,  as 
described  in  the  preceding  sections,  the  following  observations 
are  made  about  its  quality: 

(1)  The  data  has  a  high  noise  content  and  requires  filtering. 

(2)  Angular  acceleration  signals  are  kinematically  inconsis¬ 
tent  with  body  axis  rotational  velocities  and  the  Euler 
orientations.  Angular  acceleration  was  reconstructed 

by  differentiating  the  body  axis  rotational  velocities. 
This  reconstruction  was  performed  in  the  frequency 
domain  with  a  fast  Fourier  transformation  algorithm. 

(3)  For  some  of  the  flight  conditions,  the  body  axis  yaw 
rate  signal  (r)  is  inconsistent  with  the  roll  rate, 
lateral  acceleration  and  sideslip  measurements.  The 
absence  of  an  airplane  heading  signal  (¥)  made  the 
identification  of  correction  terms  for  R  difficult. 

(4)  The  angle-of-attack  measurement  requires  both  a  scale 
factor  and  bias  correction  term  to  make  the  side  force 
equation  consistent  (i.e.,  &  =  -R  +  (tanot)P  +  A  /V- 

+  g  cose  sin4)/V,p)  .  y  1 

3.3  MODEL  STRUCTURE  DETERMINATION 


System  identification  is  a  systematic  methodology  for  analyzing 
test  data  to  obtain  maximally  accurate  dynamic  equation  coefficient 
estimates.  Its  application  to  multivariable  nonlinear  systems, 
however,  requires  an  additional  step  before  actual  identification  of 
parameters  is  attempted.  This  step  is  a  mathematical  model  deter¬ 
mination  stage  where,  of  all  possible  nonlinear  aerodynamic  effects 
postulated,  those  of  greatest  significance  are  isolated  and  retained 
for  further  analysis.  There  are  two  basic  reasons  for  this  reduction 
of  order.  First,  a  complete  nonlinear  model  places  severe  computa¬ 
tional  demands  on  parameter  identification  algorithms;  hence,  the  model 
generally  must  be  reduced  in  order.  Secondly,  if  all  possible  aero¬ 
dynamic  contributions  are  retained,  the  problem  is  "overparameterized" 
by  the  inclusion  of  terms  which  are  not  identifiable.  This  may  produce 
divergence  of  the  identification  algorithm.  The  following  subsections 
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present  an  overview  of  the  model  structure  determination  procedure, 
the  formulation  of  the  T-2  aerodynamic  structure,  and  the  selected 
model  structure. 

3.3.1  Overview  of  Model  Structure  Determination 


The  model  structure  determination  procedure  is  illustrated 
in  Figure  3.3;  it  includes  the  following  tasks:  model  structure 
hypothesis,  data  processing  to  define  the  most  significant  model 
parameters,  and,  finally,  model  validation. 

3. 3. 1.1.  Model  Structure  Hypothesis 

Wind  tunnel  data,  theoretical  data  and  a  general  understanding 
of  aircraft  aerodynamic  characteristics  comprise  the  a  priori  infor¬ 
mation  on  which  functional  models  were  first  based  in  this  effort. 
While  the  aerodynamic  characteristics  of  the  T-2C  contain  no 
significant  surprises,  such  information  as  the  stall  angle  of 
attack,  the  shape  of  the  post-stall  lift  curve,  drag  level  and 
polar  shape,  and  the  basic  nonlinear  characteristics  of  stability 
and  damping  derivatives  provided  useful  guidance  in  the  initial 
formulation  of  the  aerodynamic  model  structures. 

Wind  tunnel  data,  properly  reduced  and  corrected,  are  the 
best  source  of  most  a  priori  aerodynamic  information.  However, 
the  effects  of  Reynolds  number,  model  support  interference,  and 
minor  configuration  changes  are  often  difficult  to  account  for  when 
correcting  the  data  from  wind  tunnel  to  full  scale.  The  best  source 
of  dynamic  derivative  estimates  (pitch  damping,  roll  damping,  etc.) 
is  theory,  since  these  parameters  cannot  be  accurately  determined 
in  a  wind  tunnel.  Theoretical  methods  are  thus  a  valuable  supplement 
to  the  experimental  data,  allowing  examination  of  the  sensitivity 
of  the  aerodynamic  data  to  configuration  and  flight  condition  param¬ 
eters.  To  the  extent  that  theoretical  methods  can  supplant  wind 
tunnel  testing,  their  continued  development  will  be  useful  to  model 
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Figure  3.3.  Elements  of  Model  Structure  Determination 
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structure  determination  methods.  T-2  theoretical  and  wind  tunnel 
aerodynamic  data  were  reviewed  in  Chapter  2. 

Model  structure  hypothesis  is  the  task  of  using  a  priori 
information  to  formulate  functional  relationships  between  the 
aerodynamic  force  and  moment  coefficients  and  the  independent 
variables  (e.g.,  a,  8,  P,  5r,  etc.).  This  process  is  illustrated 
by  the  following  example  for  the  yawing  moment  coefficient  equation. 


(1)  Select  the  independent  variables: 

Cn  =  f  (ct ,  8,  P,  R,  6r,  <5a) 

or  i.e.,  yawing  moment  coefficient  is  a  function  of 
angle-of-attack ,  sideslip,  roll  and  yaw  rate,  rudder 
and  aileron  position 

(2)  Partition  independent  variables  into  reasonable  groups: 

C  -  AC  +  AC  +  AC  +  AC 

nS I DESLIP  "DYNAMIC  nRUDDER  "AILERON 

(3)  Select  functional  relationships  for  each  group: 


AC. 


n 


SIDESLIP 


(K 


N1 


+  KN2a  +  KN3Aal  + 


KN4Aa2)e  +  KN5Ag 

To  effectively  handle  a  nonlinear  aerodynamic  model,  such  as 
that  required  for  the  T-2C,  a  spline  model  form  is  adopted.  The 
spline  terms  (i.e.,  Aa^ ,  Aa^,  AB  in  the  above  equation)  are  piece- 
wise  continuous  polynomials.  The  nature  of  the  spline  approach 
is  graphically  defined  in  Figure  3.4.  Because  the  spline  terms  are 
defined  over  a  restricted  region,  two  modeling  benefits  result: 

(1)  only  lower  order  terms  are  required,  and  (2)  prediction  error 
for  model  extrapolation  is  low.  Spline  models  have  the  additional 
advantage  that  they  can  be  formulated  directly  from  an  understanding 
of  aircraft  characteristics  as  shown  in  the  example  (Figure  3.4)  and 
their  calculation  (or  influence)  can  be  controlled  by  logic 
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•  A  PRIORI  DATA 
WIND  TUNNEL 
CHARACTERISTICS 

6*0 


6 


•ate  SPLINE  TERN  DEFINITIONS 
C  IF  a<a.j 
a  -  Oj  IF  a>a. 

£  +  £  IF  6<  -  6i 
<ia  *  0  IF  !6j<6i 

6  -  61  IF  6>£1 


whtre 

Kn1  IS  THE  ESTIMATE  OF  Cn  FOR  a  »  6  *  0 

IS  THE  ESTIMATE  OF  CN  FOR  a  t  0,  6  *  0 

Kjj.  CORRECTS  THE  ESTIMATE  OF  C„  FOR  ,  £  •  0 
^  CORRECTS  THE  ESTIMATE  OF  CN  FOR  a>a2  l  «  0 

Kjjg  CORRECTS  THE  ESTIMATE  OF  CN  FOR  £  t  0 


Figure  3.4.  Spline  Equation  Formulation  Example 
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variables.  For  example,  to  account  for  hysteresis,  the  inclusion 
of  a  spline  term  can  be  controlled  by  the  value  of  da/dt. 

3. 3. 1.2  Model  Structure  Data  Processing 

After  the  aerodynamic  force  and  moment  coefficients  have  been 
hypothesized,  the  next  task  is  to  determine  which  of  the  parameters 
are  significant.  This  is  done  by  reconstructing  the  aerodynamic 
force  and  moment  coefficients  from  flight  test  measurements  and 
then  using  an  equation  error  identification  algorithm  to  search  for 
the  most  significant  parameters.  The  following  equation  illustrates 
how  an  aerodynamic  moment  coefficient  (C  )  is  reconstructed  from 
flight  measurements: 

cn  ■  ^  t!z  at  -  :xz  'at  -  W  *  ']y  -  VQI 

q  =  dynamic  pressure 
Sw  *  wing  area 

bw  *  wing  span 

P,Q,R  =  roll,  pitch  and  yaw  rates 

jg-  ()  *  time  derivative  of  () 

Ix>Iy»Iz  =  r°H»  pitch,  yaw  inertias 

I  ■  product  of  inertia 

X  Id 

It  can  be  seen  from  this  equation  that  if  any  of  the  flight 
measurements  are  biased  or  have  a  scale  factor  error,  the  recon¬ 
structed  aerodynamic  coefficient  and  the  resulting  model  parameter 
estimates  will  also  be  biased. 

An  optimal  subset  regression  (OSR)  algorithm  is  used  to 
identify  the  model  structure  from  the  reconstructed  aerodynamic 
coefficient.  OSR  is  an  algorithm  which  adds  and  deletes  variables 
to  a  particular  model  in  an  iterative  manner.  The  algorithm  uses 
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statistical  hypothesis  testing  techniques  based  on  the  Fisher  F 
ratio  (e.g.,  F-tests)  .  Formally,  this  ratio  measures  the  difference 
in  fit  error  with  the  current  model  relative  to  the  error  due  to 
noise  and  model  uncertainties.  A  "total"  F-ratio  measures  the 
entire  model  fit  relative  to  error  and  a  "partial"  F-ratio  measures 
the  incremental  improvement  in  fit  due  to  addition  or  deletion  of 
a  parameter  in  the  model.  A  generalized  flow  chart  is  shown  in 
Figure  3.5.  Starting  with  a  list  of  possible  variables,  the  algorithm 
enters  the  first  variable  with  the  highest  partial  correlation  to 
the  observations  y.  The  contribution  of  this  variable  to  reducing 
the  fit  error  is  made,  and  a  new  variable  entered.  Subsequent 
tests  add  and  delete  variables  to  improve  the  "fit".  The  final 
subset  of  0  which  results  from  the  procedure  is  one  within 
confidence  bounds  set  by  the  user  (say,  95  percent  or  99  percent). 

The  model  structure  selected  from  the  OSR  results  is  based  on  a 
consideration  of  two  statistical  figures  of  merit  and  a  comparison 
of  candidate  model  structures  with  wind  tunnel  and  theoretical 
data.  This  comparison  also  provides  a  means  to  validate  the  selected 

model  structure.  The  two  statistical  figures  of  merit  are  the 

? 

multiple  correlation  coefficient  (R  )  and  the  equation  F  ratio. 

2  2 

For  a  perfect  fit,  R  =  1.  Generally,  R  -*1  as  additional  terms 

are  added  to  the  model.  The  equation  F  ratio  relates  fit  goodness 

to  fit  error  weighted  by  the  degrees  of  freedom  of  the  model.  The 

2 

desired  model  is  one  that  has  a  good  fit  (R  ->-1)  with  the 
equation  F  ratio  maximized, 

3.3.2  T-2  Aerodynamic  Formulation 


Using  the  T-2  wind  tunnel  and  theoretical  aerodynamic  data 
base  (Chapter  II)  and  the  model  hypothesis  principles  discussed  in 
Section  3. 3. 1.1,  generalized  aerodynamic  models  were  defined.  The 
model  for  each  of  the  aerodynamic  force  and  moment  equations  is 
presented  in  Figure  3.6  with  Table  3.2  containing  a  definition  of 
terms  from  these  equations.  Modeling  features  are  noted  as  follows: 
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LATERAL-DIRECTIONAL  MODELS 


Figure  3.6.  (Concluded) 


TABLE 


1. 


2. 


3.2.  AERODYNAMIC  MODEL  VARIABLES 


ANGLE  OF  ATTACK  SPLINES 
Ac*i  =  a  -  a-j  ^  0 

Aag  “  a  -  <*2  —  0 

Aa^  *  a  -  <*3  >_  0 

where  a-j  =  0.1745  rad  (10°) 

a2  =  0.2793  rad  (16°) 

ot-  =  0.3840  rad  (22°) 

SIDESLIP  SPLINE 

SB  -  B1  IF  6>B1 
0  IF  |B|<61 
B  +  61  IF  6<B1 

NORMALIZED  ROTATIONAL  RATES 
Ps=  s  *  (cosa  P  -  sina  R) 


(sina  P  +  cosa  R) 
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(1)  Nonlinearities  in  a  and  3  were  incorporated  with 
spline  terms.  Some  of  the  a  spline  terms  in  the  lift 
and  pitching  moment  equations  are  defined  separately 
for  airplane  nose  up  and  down  maneuvers. 

(2)  The  vertical  plane  force  coefficients  are  modeled  by 
lift  (Cl)  and  drag  (Cj))  instead  of  axial  (Cy) 
and  normal  (C^)  force  coefficients  since  the  curve 
shape  for  (Cp  is  much  simplier  than  that  for  C^. 

(3)  The  equations  included  cross  axis  coupling  terms.  The 
lift-drag  and  pitching  moment  models  included  sideslip 
and  rudder  deflection.  The  yawing  moment-rolling 
moment  and  sideforce  equations  are  dependent  on  angle- 
of  attack. 

The  aerodynamic  models  described  in  Figure  3.6  were  implemented 
in  the  OSR  algorithm  and  the  flight  measurements  for  the  test 
conditions  described  in  Table  2.3  were  processed.  Since  these  flight 
conditions  represented  different  kinds  of  maneuvers,  the  most 
significant  model  structure,  as  defined  by  OSR,  varied  from  run 

r 

to  run.  Table  3.3  shows  the  correlation  coefficients  obtained 

for  each  equation  and  flight  condition.  Its  apparent  that  the 

information  content  of  some  maneuvers  was  much  better  than  that 

for  others.  Also  the  models  for  C  were  worse  than  those  for  the 

n 

other  force  and  moment  equations.  The  problem  in  modeling  Cn 
is  that  the  yaw  rate  measurement,  which  is  used  to  reconstruct 
Cn,  was  a  poor  signal  (see  Section  3.1).  By  interpreting  the  OSR 
results  with  much  engineering  judgment,  the  results  were  combined 
to  produce  composite  longitudinal  and  lateral -direct ional  models. 
These  are  presented  in  Tables  3.4  and  3.5. 


3.3.3  Model  Structure  Validation 


The  selected  aerodynamic  models  were  validated  by  comparison 
to  wind  tunnel  data  and  theoretical  prediction.  Figures  3.7  to  3.12 
show  the  results  of  this  comparison. 

Figures  3.7  to  3.9  show  C L,  Cm>  and  Cp.  Of  interest  in  these 
curves,  particularly  in  C^,  is  an  aerodynamic  hysteresis  detected 


TABLE  3.3.  STATISTICAL  SIGNIFICANCE  OF  MODEL 
STRUCTURE  DETERMINATION4 


4  RAMP  *  SINE 
e 

AFT  *  PULSES 

LIMIT  CYCLE 
«e  OOUBLET 

«e  RANDOM 
«A  RANOOM 

SEQUENTIAL  00U8LETS 
6r  RANOOM 

SEQ.  RANOOM.  6^0“ 
SEQ.  RANOOM,  Bj'-S'’ 
SEQ.  RANOOM,  8t*50 
«e  AFT,  FULL  «R 
«  AFT,  FULL  6, 

C  O 

SHALLOW  BANK  TO 
STALL 

STEEP  8ANK  TO  STALL 
PULL-UP  TO  STALL 
CO-ORD.  SPIN  ENTRY 


Variance  explained  by 


CL 

CD 

.56 

.98 

.19 

.74 

.68 

.96 

.97 

.98 

.96 

.98 

.80 

.82 

.97 

.98 

.77 

.84 

.96 

.96 

.81 

.94 

.88 

.94 

.86 

.96 

.72 

.93 

.48 

.98 

.89 

.98 

.95 

.97 

.80 

.96 

AERODYNAMIC  COEFFICIENT 


.72  .06 

.58  .04 


.98  .96  .78  .18 

.82  .59  .99  .66 

.98  .94  .99  .82 

.84  .74  .96  .84 

.96  .96  .68  .19 

.94  .92  .91  .48 

.94  .94  .91  .55 


.93  .83 


.84 


.48  .04 


1  is  a  perfect  explanation) 


KVW-- 


31**^*. T,wr.'  ^ 


TABLE  3.4.  IDENTIFIED  LONGITUDINAL  MATH  MODEL  STRUCTURES 
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Figure  3.11.  Variation  of  Elevator  Effec-  Figure  3.12.  Variation  of  CR 

tiveness  with  Angle  of  Attack 

Angle  of  Attack 


by  allowing  the  coefficients  of  Aa  to  assume  different  values 

depending  on  whether  pitch  rate  was  positive  or  negative.  This 

effect  is  in  general  agreement  with  previous  studies  of  hysteresis 

but  it  is  believed  that  this  is  the  first  time  it  has  been 

identified  from  flight  data.  The  lift  deficit  causes  a  similar 

shift  in  C  over  the  same  a  range, 
m 

The  Cm  results.  Figure  3.8,  show  that  static  stability 
Cj^,  has  been  closely  identified.  An  offset  in  trim  angle  of 
attack  may  be  caused  by  an  elevator  calibration  error,  as  show 
on  the  figure.  The  pitching  moment  pitch  up  at  a  s  20°,  likely 
due  to  wing- tail  interference,  was  not  identified  due  to  the 
limited  number  of  spline  segments  available. 

Drag  correlation,  shown  in  Figure  3.9,  was  hampered  by  lack  of 
data  on  engine  thrust  level.  This  is  needed  because  the  net  X-axis 
force,  T-D,  is  identified  from  acceleration  data,  and  thrust  must 
be  known  to  determine  drag.  In  this  case,  the  thrust  was  computed 
from  estimated  drag  in  the  trim  condition,  using  wind  tunnel  data. 
More  accurate  ax  instrumentation,  engine  parameter  instrumenta¬ 
tion  (e.g.,  RPM) ,  and  models  for  installed  thrust  are  required 
in  order  to  precisely  define  the  shape  of  the  drag  polar. 

Figures  3.10  to  3.22  show  other  control  effectiveness,  aero¬ 
dynamic  damping,  and  aerodynamic  force  coefficients  and  derivatives 
identified  from  flight  data,  compared  to  values  from  other  sources. 
In  general,  identified  values  were  used  in  preference  to  other 
data,  unless  levels  of  confidence  were  very  low. 

Figures  3.14  and  3.15  show  in  detail  the  agreement  between 
the  nonlinear  models  and  wind  tunnel  data,  for  C„  and  C  as 
functions  of  8  at  various  values  of  a.  The  agreement  is,  in 
general,  quite  good. 

Figures  3.7  to  3.21  show  the  forms  in  which  starting  parameter 
values  were  input  to  the  nonlinear  maximum  likelihood  parameter 


ANGLE  0*  ATTACK.  0E6. 
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NOTES: 

•  POINTS  SH0U  WIND  TUNNEL  DATA 

•  8*0  . 

•  RN  *  6x10"/FT 

•  M-0.2 

Figure  3. 

13.  Variation  of  C.  with  Angle  of  Attack 

p 

ANGLE  OF  ATTACK,  OEG. 

CY  •  -0.010  ■ 

p 

PER 

OEGREE  -0.012 


-  POINTS  SHOW  W1N0  TUNNEL  DATA 

-  0-0 

-  RN  -  6x106/FT 


Figure  3.14.  Variation  of  CY  with  Angle  of  Attack 


'  ;t,0S*  1*1 


Figure  3.18.  Variation  of  C1  with  Angle  of  Attack 

r 
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mu  v  attack,  du. 


Figure  3.19.  Variation  of  Rudder  Effectiveness 
with  Angle  of  Attack 


AMU  OF  ATTACK,  DU. 


Figure  3.20.  Variation  of  Aileron  Effectiveness 
with  Angle  of  Attack 
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identification  program.  Specific  spline  regions  were  called  out 
for  which  this  program  would  attempt  to  obtain  improved  parameter 
values,  in  the  sense  of  improved  overall  model  fit.  The  results 
of  this  investigation  are  discussed  in  the  following  section. 

3.4  PARAMETER  IDENTIFICATION  AND  VERIFICATION 
3.4.1  Approach 

Two  operations  are  now  performed  with  the  aerodynamic  coeffi¬ 
cients  determined  above:  (1)  the  coefficients  are  used  as  starting 
values  in  a  maximum  likelihood  parameter  estimation  analysis 
searching  for  improved  estimate  values;  and  (2)  dynamic  time  his¬ 
tories  are  generated  from  the  model  equations  and  compared  to  the 
measurements,  giving  an  indication  of  model  fidelity. 

The  6  degree  of  freedom  model  is  divided  into  the  customary 
3  DOF  longitudinal  and  lateral-directional  models  in  order  to 
simplify  error  interpretation  and  reduce  the  number  of  parameters 
to  be  identified  at  one  time.  Lateral -directional  variables  are 
supplied  by  lateral -directional  measurements  read  in  during 
the  solution  of  the  longitudinal  equations  of  motion.  This  does 
not  prohibit  the  estimation  of  coupling  derivatives,  such  as  the 
dependence  of  a  longitudinal  parameter  on  a  lateral-directional 
degree  of  freedom,  C^g,  f°r  example>  since  8  is  read  in  to 
the  longitudinal  equations  in  the  same  manner  as  a  control  deflection. 

Several  runs  were  selected  for  further  analysis,  based  on  the 
range  of  a  and  8  available  and  on  the  general  quality  of  the 

correlation  coefficient  characteristics  established  during 
the  model  structure  determination  phase  of  the  study  (see  Table 
3.3).  The  runs  selected  are,  for  the  longitudinal  case,  the  random 
elevator  and  elevator  doublet  maneuvers,  and  for  the  lateral- 
directional  case  the  random  aileron  maneuver.  The  former  longi¬ 
tudinal  maneuver  was  used  for  identification,  the  latter  for  a 
prediction  test,  the  "acid  test"  of  a  model;  the  lateral-directional 
maneuver  was  used  for  identification  only,  due  to  generally  poorer 
measurement  data  quality  (see  Section  3.2.4). 
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The  criteria  employed  in  model  fidelity  analysis  are  both 
qualitative  and  quantitative.  Obviously,  a  perfect  time  history 
match  on  a  prediction  run  is  indicative  of  a  good  model.  Recall 
that  a  prediction  run  is  one  in  which  the  control  inputs  from  one 
maneuver,  say  B,  are  applied  to  the  math  model  derived  from 
another  maneuver,  say  A,  and  the  resulting  model  motion  calcu¬ 
lations  compared  to  the  actual  cehavior  of  B.  If  perfect  agree¬ 
ment  is  seen,  the  model  has  been  shown  to  be  general  and  not  the 
result  of  a  unique  convergence  to  data  from  a  particular  maneuver. 

The  identification  process  yields  calculations  of  both  noise 
covariance,  or  dispersion  about  the  mean,  and  model  error  covari¬ 
ance,  or  the  difference  between  the  model  predictions  and  the 
measured  data.  Statistical  results  showing  large  error  covariances 
are  often  very  apparent  in  poor  time  history  matches  of  a  given 
degree  of  freedom. 

The  results  in  the  next  two  subsections  illustrate  the  type 
of  treatment  employed  in  this  study,  and  the  results  obtained. 

3.4.2  Longitudinal  Results 

Model  parameter  values  obtained  from  the  model  structure 
determination  process  were  used  as  starting  values  to  the  maximum 
likelihood  parameter  estimation  program  SCIDNT  and  improved 
estimation  of  selected  parameters  was  sought.  The  level  of 
improvement  was  evaluated  by  comparing  model  error  covariance 
magnitudes  before  and  after  the  operation,  and  by  time  history 
comparisons  of  the  final  model  to  the  measured  data. 

The  initial  longitudinal  model  parameters  were  shown  in 
Table  3.4.  Data  from  the  random  elevator  maneuver  (see  Figures 
3.22-3.26)  were  processed  for  model  identification.  The  improve¬ 
ment  in  the  model  resulting  from  running  the  maximum  likelihood 
program  is  seen  by  comparing  the  initial  model  error  covariances 
with  the  final  error  covariances,  as  shown  in  Table  3.6.  These 
values,  computed  from  time  history  comparisons  of  model  estimates 
of  the  measurements  versus  actual  data,  sho\</  that  the  model  fit 
is  much  improved  after  the  adjustment  of  parameter  values. 
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TABLE  3.6.  MEASUREMENT  MEAN  SQUARE  ESTIMATION  ERROR: 
RANDOM  ELEVATOR  MANEUVER 


a  ~ 

rad 

V  ~ 

fps 

INITIAL 

a 

V  Q 

e 

Q  ~ 

rad/sec 

\  e  - 

rad 

a 

.0255 

0  0 

0 

V 

0 

97.6  0 

0 

V  WITH 

PARAMETER  VALUES 

Q 

0 

0  .0085 

0 

FROM  MODEL  STRUCTURES 

V  DETERMINATION  ANALYSIS 

e 

0 

0  0 

.0625_ 

/ 

FINAL 

\ 

a 

.00323 

0  0 

0 

V 

0 

60.75  0 

0 

>  WITH  MAX.  LIKELIHOOD  IMPROVEMENT  ! 

Q 

0 

0  .0065 

0 

i 

0 

e 

0 

0  0 

.00072 

/  a 

.0568  rad 

'  V 

7.79  fps 

Q 

.0806  rad/ sec 

6 

.027  rad 
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The  angle-of-attack  time  history,  Figure  3.23,  shows  good  agreement 
on  frequency  and  amplitude  over  the  duration  of  the  maneuver. 
Deviations  near  peak  amplitude  contribute  the  most  toward  rms 
error  (cr)  ,  and  small  differences  in  phase  appear  as  large  differences 
in  magnitude.  Therefore,  the  overall  fit  appears  better  than  the 
oa  =  3.26°  value  might  imply. 

The  agreement  in  pitch  rate.  Figure  3.24,  an  important  pilot 
cue  and  handling  qualities  determinant,  is  excellent.  Elevator 
effectiveness  and  short-period-related  aerodynamic  parameters 
appear  to  have  been  well  identified.  Again  the  Og  =  4.62°  /sec 
value  is  misleading  relative  to  the  adequacy  of  the  fit. 

It  is  concluded  here  that  simple  measures  of  fit  accuracy, 
such  as  standard  deviation,  are  unreliable  for  judging  models,  and 
must  be  accompanied  by  engineering  judgment  and  more  comprehensive 
statistical  criteria,  such  as  correlation. 

Pitch  attitude,  Figure  3.25,  also  shows  excellent  agreement 
throughout  the  maneuver.  The  agreement  found  with  velocity, 

Figure  3.26,  is  not  good,  and  is  in  fact  worse  than  ov  *  7.8  fps 
might  imply.  The  generally  poor  quality  of  longitudinal  accelero¬ 
meter  data  is  held  responsible  for  this  lack  of  agreement. 

The  conclusion  is  that  the  model  determined  in  general  form 
by  regression  analysis  of  numerous  maneuvers,  then  refined  by 
the  parameter  identification  program,  predicts  aircraft  behavior 
very  well  in  the  nonlinear,  post-stall  flight  regime,  as  evidenced 
by  the  results  of  this  first  maneuver  analysis.  To  test  model 
generality,  a  second  maneuver  was  selected,  also  with  a  wide 
angle-of-attack  range,  as  shown  in  its  histogram,  Figure  3.27, 
but  more  in  the  linear  below-stall  regime.  The  same  model,  with 
the  control  inputs  of  the  new  maneuver,  produced  the  error 
statistics  shown  in  Table  3.7  and  the  time  history  comparisons 
shown  in  Figures  3.28  through  3.31. 

Again,  time  history  comparisons  are  found  to  be  very  good. 
Angle-of-attack  and  pitch  rate  show  well-modeled  short  period 


i#  DOUBLET  MANEUVER 

NUMBERS  ON  the  PLOT  SHOW  THE  NUMBER  OF  DATA 
"01 NTS  FOR  INDICATED  /•  REGION 


Figure  3.27, 


a/S  Histogram  for  First  Longitudinal 
Model  Prediction  Maneuver 


TABLE  3.7.  MEASUREMENT  MEAN  SQUARE  ESTIMATION  ERROR 
ELEVATOR  DOUBLET  MANEUVER 


Figure  3.28.  Prediction  of  5  Doublet  Maneuver  with 
Sp  Random  Model?  Angle  of  Attack 


Figure  3.29.  Prediction  of  6  Doublet  Maneuver 
with  6^  Random  Model:  Pitch  Rate 


Figure  3.30.  Prediction  of  6g  Doublet  Maneuver 

with  Random  Model:  Pitch  Attitude 


Figure  3.31.  Prediction  of  6  Doublet  Maneuver 
with  <$e  Random  Model:  Airspeed 
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and  control  response  characteristics;  pitch  attitude,  accordingly, 
is  good;  and  velocity,  still  the  poorest,  shows  the  proper  trend 
and  is  close  on  response  frequency. 

To  further  test  the  validity  of  this  longitudinal  model,  a 
maneuver  with  large  excursions  in  both  angle  of  attack  and  side¬ 
slip  was  selected  for  comparison.  Before  computing  time  histories, 
the  sideslip  information  was  processed  in  an  identification  mode 
to  refine  estimates  of  8-dependent  parameters.  A  histogram  of 
this  maneuver,  showing  a  and  8  excursions,  is  shown  in  Figure 
3.32.  Figures  3.33  through  3.36  show  the  time  history  comparisons 
for  this  case.  Again,  though  peak  amplitudes  are  underpredicted, 
fair  agreement  between  model  and  measurement  is  seen. 

The  conclusion  reached  from  this  longitudinal  analysis  is 
that  a  general  model  has  been  developed  in  the  approach  described, 
that  gives  results  with  accuracy  adequate  for  simulating  piloted 
handling  qualities  in  the  nonlinear,  high-a  regime. 

3.4.3  Lateral-Directional  Results 


The  lateral-directional  parameters  identified  by  regres¬ 
sion  analysis  were  shown  in  Table  3.5.  Improved  estimates  of 
these  parameters  were  sought  using  maximum  likelihood  analysis. 

The  parameters  selected  for  identification  are  shown  in 
Table  3.8.  These  parameters  are  those  judged  most  important  to 
lateral-directional  dynamics.  Table  3.9  shows  model  error  covar¬ 
iances  and  standard  deviations  before  and  after  the  maximum  likeli- 
hodd  identification  operation.  While  the  identified  parameters 
do  not  show  large  changes  in  value  from  the  starting  values,  and 
none  changed  sign,  the  error  covariances  are  reduced  for  all  but 
the  roll  rate  degree  of  freedom.  Additional  analyses  are  indicated 
to  further  reduce  the  error,  beginning  with  the  addition  of  more 
variables,  particularly  Cn  ,  to  the  identified  set. 
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NOTES: 

1)  SEQUENTIAL  DOUBLET  (6  ,  6  f  fiR)  MANEUVER 

2)  NUMBERS  ON  THE  PLOT  SHOW  THE  NUMBER  OF 
DATA  POINTS  FOR  INDICATED  a/B  REGION 
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Figure  3.32.  a/6  Histogram  for  2nd  Longitudinal 
Model  Prediction  Maneuver 
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TABLE  3.8.  COMPARISON  OF  STARTING  AND  IDENTIFIED  LATERAL 
PIRFCTIONAL  PARAMETERS 


PARAMETER* 

STARTING 

FINAL 

% 

-.0790 

-.0793 

% 

-.502 

-.598 

C. 

-.690 

-.715 

pa 

^UpAal 

5.85 

6.39 

CV 

.292 

.268 

C£ 

-.212 

-.201 

V 

.086 

-.122 

.122 

-.112 

S 

r 

-.0526 

-.018 

CY 

B 

-.338 

-.334 

*  Other  lateral-directional  para¬ 
meters  are  as  defined  in  Table  3.5 


TABLE  3.9.  MEASUREMENT  MEAN  SQUARE  ESTIMATION  ERROR: 
RANDOM  AILERON  MANEUVER 


INITIAL  MODEL  ERROR  COVARIANCE 

a 

P 

To 354  0  0  0  “ 

P 

.188  rad/sec 

r 

0  .0188  0  0 

r 

.137  rad/sec 

0  0  2.854  0 

1 .689  rad 

6 

000  .00413 

6 

.064  rad 

FINAL  MODEL 

ERROR  COVARIANCE 

a 

P 

To641 

0 

0 

0  ~ 

P 

.25  rad/sec 

r 

0 

.00252 

0 

0 

r 

.050  rad/sec 

<t> 

0 

0 

.0490 

0 

.221  rad 

8 

0 

0 

0 

,00274 

6 

.C52  rad 
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The  maneuver  selected  for  analysis  here  is  illustrated  by 
the  histogram  in  Figure  3.37.  It  represents  the  response  to  a 
random  aileron  control  input  at  stall.  Time  history  model- 
measurement  comparisons  are  shown  in  Figures  3.38  through  3.41, 
for  the  variables  p,  r,  <p ,  and  8.  In  general,  the  agreement 
is  reasonable.  Frequencies  are  well  represented,  although  some 
underprediction  of  peak  response  amplitudes  is  generally  seen. 
As  the  model  stands,  the  basic  lateral-directional  dynamics  are 
well  represented. 

NOTES: 

1 )  6a  RANDOM  MANEUVER 

2)  NUMBERS  ON  THE  PLOT  SHOW  THE  NUMBER  OF  DATA 
POINTS  FOR  INDICATED  a/8  REGION 
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Fieure  3.37. 


a/8  Histogram  for  Lateral/Directional 
Model  Identification 
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Lateral-Directional  Model  Identification 
6  Random  Maneuver:  Roll  Attitude 
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Figure  3.41.  Lateral-nirectiona 1  Model  Identification 
<5  Random  Maneuver:  Sideslip 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


4.1  CONCLUSIONS 


This  study  represents  one  of  the  first  thorough  applications 
of  advanced  nonlinear  system  identification  techniques  to  the 
problem  of  determining  aircraft  aerodynamic  characteristics  in 
the  maximum  lift  flight  regime,  and  has  produced  results  con¬ 
firming  the  capability  of  these  methods  to  develop  realistic  and 
useful  dynamics  models  for  flight  test  data  and  to  provide  a 
valuable  guide  for  data  acquisition  system  specification  in  future 
flight  test  programs .  This  study  also  demonstrated  the  value  of 
employing  wind  tunnel  data  and  theoretical  (or  empirical)  aero¬ 
dynamic  analyses  during  the  model  structure  determination  stage 
in  order  to  begin  with  a  basic  familiarity  with  the  nature  of 
nonlinearities  to  be  encountered. 

The  following  are  the  principal  conclusions  drawn  from  this 

work: 

(1)  Useful  and  accurate  system  identification  work  --  that 
is,  model  structure  determination  and  parameter  esti¬ 
mation  --  can  be  performed  even  with  poor  quality  data, 
due  to  the  ability  to  reconstruct  lost  data  and  filter 
data  with  high  noise  content.  Data  reconstruction,  by 
computation,  differentiation  and  integration,  stands 
out  as  an  almost  indispensable  tool  for  the  initial 
processing  of  flight  data.  Examination  of  the  kine¬ 
matic  consistency  of  the  data  gives  early  indication 

of  data  acquisition  errors. 

(2)  The  most  effective  approach  to  model  structure  deter¬ 
mination  employs  semi-empirical  aerodynamic  methods 
(e.g.,  DATCOM) ,  available  wind  tunnel  data,  and  aero¬ 
nautical  experience  to  establish  a  basic  model  hypo¬ 
thesis  with  the  minimum  number  of  independent  variables 
and  the  best  estimate  of  the  natures  of  the  nonlinear¬ 
ities  expected  in  the  data. 


(5)  The  multiple  regression  determination  of  model  struc¬ 
tures  was  highly  effective  in  this  effort,  and 
identified  coefficients  of  significance  with  good 
numerical  accuracy,  as  judged  by  dynamic  simulation 
results. 

(4)  Separation  of  aircraft  equations  of  motion  into 
independent  longitudinal  and  lateral -direct ional  sets 
and  identifying  the  parameters  of  each  set  separately 
was  found  to  produce  good  identification  of  coeffi¬ 
cient  values  while  decreasing  computational  complexity 
and  computer  run  time. 

(5)  The  generality  of  the  models  developed  in  this  study 
was  demonstrated  by  predicting  the  time  history  of 
one  control  response  maneuver  using  a  model  developed 
from  a  different  maneuver.  This  generality  resulted 
from  the  substantial  effort  taken  to  define  model 
structures  containing  aerodynamically  relevant  non¬ 
linear  terms,  basing  them  on  the  combined  results  of 
the  analysis  of  many  flight  maneuvers  and  on  experience 

(6)  The  aerodynamic  coefficients  identified  in  this  effort 
agreed  well  in  value  with  wind  tunnel  data  and  theoret¬ 
ical  aerodynamic  calculations  of  derivatives  in  the 
post-stall  angle  of  attack  regime.  Iterative  appli¬ 
cations  of  the  model  structure  determination  procedure 
enabled  the  optimization  of  spline  function  break  point 
to  the  nonlinearities  exhibited  by  the  coefficients 
within  the  range  of  angle  of  attack  covered  by  the 
data.  Since  the  majority  of  the  data  were  in  the 
post-stall  regime,  identification  of  coefficients 

in  the  low-a,  linear  regime  were  expectedly  of  lower 
accuracy. 

(7)  The  numerical  approach  to  data  processing  in  this 
system  identification  methodology  resulted  in  the 
identification  of  such  highly  nonlinear  aerodynamic 
characteristics  as  hysteresis  in  lift,  drag  and 
pitching  moment  around  the  stall  angle  of  attack, 
and  dynamic  deviations  with  sign  changes  and  reflex 
curvature  versus  angle  of  attack. 

(8)  The  aerodynamic  coefficients  extracted  from  the  given 
flight  data  comprise  a  basic  data  set  for  a  flight 
simulator  math  model  and  would  be  adequate  to  conduct 
piloted  simulations  of  vehicle  flying  qualities  in  the 
high  angle  of  attack  regime. 
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The  results  of  this  effort  are  felt  to  be  highly  significant 
in  that  a  unified  procedure  has  been  developed  which  employs  all 
of  flight,  wind  tunnel,  and  theoretical  techniques  to  achieve 
improved  understanding  of  aircraft  characteristics.  The  ability 
to  develop  accurate  models  of  aircraft  behavior  from  flight 
measurements  presents  the  opportunity  to  plan  an  aircraft  program's 
wind  tunnel,  theoretical,  and  flight  test  work  in  a  coordinated 
fashion  so  that  each  activity  fills  its  role  with  maximum  effec¬ 
tiveness,  and  the  best  capabilities  of  each  are  used  to  their 
full  and  timely  advantage.  The  success  of  the  system  identifi¬ 
cation  procedure  demonstrated  in  this  effort  signifies  that 
flight  data  analysis  may  now  be  performed  in  as  rigorous  and 
useful  a  fashion  as  other  experimental  and  theoretical  methods 
have  been  used  in  the  past.  As  a  guide  to  theoretical  methods 
development  alone,  through  its  ability  to  verify  predictions  of 
the  dynamic  stability  derivatives  of  increasing  importance  in 
advanced  aircraft  mission  performance,  the  system  identification 
method  can  furnish  information  of  value  to  a  wide  range  of  air¬ 
craft  development  programs.  And,  as  future  flight  test  appli¬ 
cations  extend  into  the  transonic  flight  regime,  the  system 
identification  method  will,  for  the  first  time,  enable  the  study 
of  the  complex  aerodynamic  phenomena  that  affect  aircraft  char¬ 
acteristics  in  that  regime  and  enhance  the  work  that  previously 
could  only  be  done  in  wind  tunnels  using  extensive  subjective 
engineering  judgment.  In  all,  the  conclusions  reached  in  this 
effort  are  encouraging  that  an  advance  in  the  state  of  the  art 
of  aircraft  engineering  has  been  achieved. 

4.2  RECOMMENDATIONS 

This  effort  has  been  a  valuable  exercise  in  the  identifi¬ 
cation  of  nonlinear  aerodynamic  characteristics,  because  it  has 
demonstrated  the  basic  viability  of  the  procedure.  The  following 


75 


recommendations  refer  to  improvements  which  would  lead  to  more 
effective  uses  of  the  system  identification  procedure  and 
enhanced  data  processing  effectiveness: 

(1)  Data  consistency  analysis  should  be  made  an  integral 
part  of  the  aircraft  parameter  identification  pro¬ 
cedure,  to  accelerate  the  process  of  data  reconstruc¬ 
tion. 

(2)  Aircraft  instrumentation  should  be  carefully  coordinated 
in  terms  of  redundancy,  response,  range  and  resolution 
to  provide  optimum  data  for  accurate  parameter  identifi¬ 
cation.  Careful  experimental  procedures  should  be  used 
to  assure  data  acquisition  system  integrity  prior  to 
flight . 

(3)  Program  development  efforts  should  be  directed  toward 
automating  the  powerful  spline-function  representation 
of  coefficients  within  the  regression  program  used  in 
model  structure  determination,  as  the  existence  of 

an  accurate  model  greatly  improves  the  accuracy  of 
later  parameter  identification  results. 

(4)  Further  work  should  be  directed  toward  developing  a 
general  nonlinear  aerodynamic  identification  model 
suitable  for  fixed-wing  aircraft  analysis.  Following 
the  approach  taken  in  this  effort,  the  model  should 
be  in  a  form  readily  implementable  in  a  flight  simu¬ 
lator  math  model. 

(5)  Subsequent  efforts  should  be  directed  toward  the 
integration  of  system  identification  (as  applied  to 
flight  test),  wind  tunnel  experimentation,  and 
theoretical  methods  into  procedures  and  tools  by 
which  they  may  most  effectively  contribute  to  the 
needs  of  aircraft  development  and  evaluation  programs. 
Commonalities  in  requirements  and  capabilities,  and 
areas  in  which  strengths  can  be  combined,  should  be 
determined.  A  research  program  in  this  area  is  almost 
certain  to  produce  interesting  and  useful  results  lead¬ 
ing  toward  a  unified  aerodynamic  treatment  of  aircraft. 

(6)  Extensions  of  this  method  to  the  determination  of 
structural  parameters  from  flight  test  data,  and  the 
study  of  their  correlation  to  aerodynamic  results, 
particularly  in  the  transonic  regime,  should  be 
pursued . 

These  recommendations,  it  is  felt,  will  eanble  the  result 
of  this  study  to  be  fully  utilized  in  improving  the  ultimate 
effectiveness  of  aircraft  engineering  activities. 
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APPENDIX  A 
T-2C  SPECIFICATIONS 


n - — ^ 

WING 

i 

Total  area  (includes  flap,  aileron  and 

39.39  ft2  covered  by  fuselage) 

254.86  ft2 

1  *» 

Net  surface  area  (vetted) 

424.85  ft2 

y 

\ A 

Span  (perpendicular  to  plana  of 
symmetry)  Including  tiptanka 

38.13  ft 

t:  ARw 

Aspect  Ratio 

5.07 

I; 

Taper  Ratio 

.496 

t- 

r  r 

■  w 

Dihedral  Angle 

+3* 

r 

r 

Chord  (in  streamline  direction) 

P  cr 

Root  (Wing  Sta.  0) 

114.20  in 

r  ct 

Tip  Chord  (Wing  Sta.  214.242) 

56.63  in 

fj 

['i 

(Equivalent) 

Mean  aerodynamic  chord 

88.88  in 

Lj 

(Wing  Sta.  95.078) 

| 

Location  of  25X  MAC 

F.S.  219.697 

V 

Sweep back  of  25X  element 

2*17' 

** 

Incidence  angle 

Root  Chord 

2* 

Tip  Chord 

-1* 

R 

t i 

Airfoil  Section  (root  and  tip  in 
streamline  direction) 

*MA  Modified 

NASA641A212 

2  -  .8*  (MOD) 

(flaps  and 
ailerons  rigged 

3*  up)  » 

| 

°v 

Race  of  Taper 

0.2671 

i 

■; 
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* 

FLAP  (Data  for  One) 

Type 

Sj  Area 

bf  Span  (perpendicular  to  plane  of  syaaetry) 

ci  Inboard  chord  (Wing  Sta.  27.09) 

c  Outboard  chord  (Wing  Sta.  127.54) 

o 

c£/c^  Ratio  flap  chord  to  vlng  chord  (erg.) 
b^/bv  Ratio  flap  span  to  ving  eeai-epan 

I 

*.  Flap  deflection,  aaxiaua  (froa  uprigged 

*  position) 

Flap  in  neutral  position 

AILERON 

Type 

S^  Area  (aft  of  hinge  line  and  including  cab) 

b^  Span  (perpendicular  to  plane  of  syaaetry) 

ct  Inboard  chord  (Wing  Sta.  128.69) 

c  Outboard  chord  (Wing  Sea.  208.26) 

o 

c  /c^  Ratio  aileron  chord  (aft  B.L.)  to  wind  chord 
b^/bv  Ratio  aileron  span  to  wind  seai-apan 

I 

8  Aileron  deflection,  maximal  (froa  neutral 

position) 

Aileron  in  neutral  position 
Aerodynamic  Balance 

S.  Balance  area  forward  of  the  H.L. 

(including  50X  of  fabric  seal) 

c./c  Ratio  balance  chord  to  aileron  chord 

D  ft 


Single  Slotted 
22.78  ft2 
101.75  in 
39.39  in 
29.63  in 
.37 
.475 

33* 

3*  Op 

Straight  Sided 
9.5  ft2 
79.57  in 
20  in 
14.66  in 
.25 
.374 

-12*  Op,  +13*  Dn 
3*  Op 

Sealed  paddle 
balance 

4.45  ft2 
.42 
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AILERON  -  (Coat'd) 

Static  balance 


Irreversible  full  power  syscea 
AILERON  TRIM  TAB 

Ground  adjustable  fixed  tab  on  each 
aileron 

Area  (each) 

HORIZONTAL  TAIL* 

S.  Total  area  (includes  3.07  ft2  covered 

by  vertical  tail  and  fairing) 

Net  area 

Net  surface  area  (vetted) 

Span 

Aspect  Ratio 
Taper  Ratio 
Dihedral  Angle 
Sveepback  of  2SI  element 
Chord  (in  streamline  direction) 

Root  (H.T.  Sta.  0) 

Equivalent  tip  chord 
(H.T.  Sta.  106.448) 

Mean  aerodynamic  chord 
(H.T.  Sta.  47.78) 

Incidence  angle 

Airfoil  section  (root  and  tip  in  streamline 
direction) 


n€th 

s 

“h 

ARh 

xh 

rh 

*h 


Weighted  paddle 
balance 

Hydraulic 


.07  ft2 


72.29  ft2 
69.22  ft2 

146.38  ft2 
17.91  ft 
4.42 
0.50 
0* 

15* 

64.61  in 
33.05  in 

50.447  in 

0* 

NASA  65A012 


★percent  lines  base  on  horizontal  prior  to  addition  of  trailing  edga  extension. 
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202.58  la 


HORIZONTAL  TAIL  -  (Coat'd) 


Tail  length  (.25  c  to  .25  c.  ) 
v  n 


HORIZONTAL  STABILIZER 


ELEVATOR 

S_ 


o 

ce/ch 


be/bh 


Area  atabilizer,  total 
Stabilizer  incidence  angle 


Total  area  (excluding  balance  area  forward 
of  the  hinge  line) 

Span  (between  equivalent  chords) 

(one  elevator  only) 

Inboard  chord  (B.P.  3.906) 

Outboard  chord  (B.P.  105.877) 

Ratio  elevator  chord  (aft  H.L.)  to 
horizontal  tail  chord 

Ratio  elevator  span  to  horizontal  tail 
span 

Elevator  deflection  maximum 
Boost:  Push  force  2.95:1 

Pull  force  2.95:1  to  8  lbs 
then  6.0:1 
Static  balance 


Aerodynamic  balance 

Sb  Balance  area  forward  of  hinge  line 

c./c  Ratio  balance  chord  to  elevator  chord 

D  • 

Nose  factor 

Point  of  taagency  for  nose  factor  is 
at  elevator  hinge  line 


42.5  ft2 
0* 


21.00  ft2 

101.97  in 
18.85  in 
10.52  in 

.310 

0.936 

27*  Op,  15*  On 
Hydraulic 


Weighted  Leading 
Edge 

Overhang 

5.72  ft2 

0.322 


8  4 


0.60 


ELEVATOR  TRIM  TAB 


S£  Aru  (each) 

b£  Span,  Equivalent  (B.P.  8.93  to  54.53) 

c£  Chord ,  constant 

bc/t>c  Ratio  tab  span  to  alevator  span 

6c  Tab  deflection 

VERTICAL  TAIL 

S  Total  area  (includes  4.38  ft2  blanketed 

9  by  fuselage  plus  2.14  ft2  blanketed  by 

horizontal  tall) 

S  _  Net  area 

net 

v 

A^  Net  surface  area  (vetted) 

Ad  Net  surface  area  of  dorsal  fin  (vetted) 

by  Span,  unblanketed 

ARy  Aspect  Ratio 

Taper  Ratio 

c?  Chord  (in  streamline  direction) 

Root  (W.P.  +  33.000) 

c£  Equivalent  Tip  Chord  (W.P.  +  129.41) 

cv  Mean  aerodynamic  chord  (W.P.  +  73.92) 

Ay  Sveepback  (25Z  chord) 

Airfoil  Section 

ly  Tail  length  (.25  cy  to  .25  c^ 

VERTICAL  FIN 

S,  Area  (including  2.14  ft2  blanketed  by 

horizontal  tail  and  excluding  dorsal  fin) 


2.36  ft2 

46.10  in 

6.5  in 

0.462  in 

L.H.  10*  Up, 
13* 

R.H.  0* 

13* 


40.33  ft2 
33.86  ft2 

79.18  ft2 
18.12  ft2 
8.04  ft 
1.80 
.375 

78.14  in 
29.38  in 
58.47  in 
30* 

NASA  63A012 
194.05  in 


29.87 


8  5 


v-Ti*- 


■JWEL 


VERTICAL  FIN  (Coat'd) 

if  Angle  with  respect  to  airplane  plane  of 

syenetry 

RUDDER 

Sr  Total  area 

Sr  Upper  surface 
u 

S  Lover  surface 
rl 

b£  Span,  equivalent 

br  Upper  surface 
u 

b  Lover  surface 
rl 

cr  Upper  chord  (W.P.  96.00) 

u 

c  Lower  chord  (W.P.  +  9.91) 

rl 

c  /c  Ratio  rudder  chord  (aft  H.L.)  to  vertical 

tail  chord 

cr/c^  upper  surface  <?  W.P.  96.00 

c/c  Lover  surface 
r  v 

*r  Rudder  deflection,  maxima 

Boost 

Aerodynamic  balance 

Sb  Balance  area  forward  of  hinge  line 

cb/cf  Ratio  balance  chord  to  rudder  chord 

cb/c  Upper  surface  @  W.P.  96.00 
u 

c, /c  Lower  surface 
d  r 

Static  balance 
Rose  factor 

Point  of  tangency  for  nose  factor  is  at 
rudder  hinge  line 


0* 

9.13  ft2 
3.23  ft2 

5.90  ft2 

31.94  in 
42.99 
12.59  in 
22.45  in 

.266 

.250 

25*  Rt.,  25*  Lt. 
Rone 

Overhang 
2.41  ft2 

.234 

.24 

Weighted  leading 

edge 

0.40 
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RUDDER  TRIM  TAB 


V»r 


FUSELAGE 

lx 


f 

L/D 

CANOPY 

lc 

Fc 

*c 

L/Dc 

NACELLE 


L/D. 


Area 

Span,  equivalent  (W.P.  14.94  to  W.P.  53.00) 

Chord,  constant 

Ratio  tab  span  to  rudder  span 

Tab  deflection,  aaximm 

Length  (actual) 

Maxima  frontal  area  (basic  fuselage) 
Maximum  width  (basic  fuselage)  F.S.  169 
Max  1mm  depth 

Basic  fuselage  over  canopy  (F.S.  169) 
Including  ducts  (F.S.  214) 

Net  surface  area 
Fineness  ratio  (actual) 

Length  (actual) 

Maxinm  frontal  area 
Net  surface  area 
Fineness  ratio  (actual) 

Length  (actual) 

Maxima  frontal  area 

Net  surface  area 

Inlet  area  (includes  gutters) 

Fineness  ratio  (actual) 


1.60  ft2 
38.06  in 
6.0  in 
.508 

7*  Rt.,  7*  Lt. 

34.58  ft 

15.75  ft2 
54  in 

88.1  in 

73.9  in 
221.11  ft2 
5.91 

19.75  ft 
3.70  ft2 

73.10  ft2 

8.8 


23.71  ft 
10.50  ft2 
206.0  ft2 
3.1  ft2 


5.025 


SPEED  BRAKE  (Data  for  one  tide  only) 


Type 

One  Piece 

Location 

Side  of  Aft 

Fuselage 

Nuber 

Two 

SJ 

Area  (Planform) 

8.00  ft2 

FJ 

Area  (frontal) 

4.24  ft2 

*j 

Maximum  deflection 

32* 

TIP  TANK 

(Data  for  one  tank  only) 

ltt 

Overall  length 

142.75  in 

dtt 

Maximum  diameter  (Tank  Sta.  61.875) 

20.00  in 

L/D 

Fineness  ratio 

7.14 

Sstt 

Side  area  (projected) 

14.1  ft2 

Sptt 

Planform  area  (projected) 

14.2  ft2 

Volume 

15.3  ft3 

\t 

Total  Surface  Area 

44.30  ft2 

A  , 
net.. 

Net  Surface  Area  (vetted) 

42.40  ft2 
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